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ABSTRACT 


The  computer  code  described  herein  permits  the  solution  of  the 
Rankine-Hugoniot  equations  for  normal  shock  waves  in  ideal  gas 
mixtures  for  a  large  number  of  optional  sets  of  conditions.  The 
code  has  been  written  partly  in  FORTRAN  II  and  partly  in  FAP,  and  is 
intended  for  use  on  the  IBM  7090  or  7094  computers  operating  under 
the  MONITOR  system. 
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INTRODUCTION 


For  many  practical  purposes,  a  shock  wave  in  a  gas  can  be  con¬ 
sidered  as  a  discontinuous  jump  between  initial  and  final  conditions, 
and  use  can  be  made  of  the  relationships  between  these  two  sets  of 
conditions.  These  relationships  are  given  for  normal  shock  waves  by 
the  so-called  Rankine-Hugoniot  equations,  coupled  with  the  thermal 
and  caloric  equations  of  state.  These  are  a  set  of  algebraic 
equations  which  in  general  require  numerical  solution  because  of  the 
complex  dependence  of  the  equations  of  state  upon  internal  molecular 
excitation  and  chemical  reaction. 

The  computer  code  which  this  report  describes  permits  the 
solution  of  the  Rankine-Hugoniot  equations  for  a  large  range  of  cases. 
The  central  problem  is  that  of  a  shock  wave  of  arbitrary  strength  in 
a  gas  at  rest.  In  addition,  the  code  permits  the  consideration  of 
the  reflection  of  such  a  shock  wave  off  a  rigid  wall  or  the  one¬ 
dimensional  detonation  of  the  initial  gas  according  to  the  Chapman- 
Jouguet  model.  The  composition  of  the  initial  gas  may  be  specified 
arbitrarily,  and  chemical  changes  occurring  between  the  initial  and 
final  states  may  be  either  prohibited,  specified  arbitrarily,  or 
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* 

determined  by  conditions  of  chemical  equilibrium  in  the  final  state* 

Provision  is  also  made  for  the  consideration  of  chemical  equilibrium 
problems  independent  of  the  occurrence  of  shock  waves. 

The  code  is  intended  for  use  on  the  IBM-7090  or  7094  computers 
under  the  MONITOR  system.  It  is  written  partly  in  FORTRAN  II  and 
partly  in  FAP. 
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Chapter  1 
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SHOCK  WAVES  IN  REACTIVE  GASES 

1.1  NORMAL  SHOCK  WAVES 

The  three  Rankine-Hugoniot  equations  express  the  conservation 
"*  of  energy,  mass,  and  momentum  in  the  flow  through  a  normal  shock 

wave.  The  ordinary  forms  of  these  equations  describe  steady,  one¬ 
dimensional,  adiabatic  flow  with  no  external  forces.  In  a  coordinate 
system  fixed  in  the  shock,  they  may  be  conveniently  written  to  state 
that  the  specific  stagnation  enthalpy  and  the  mass  flux  are  the  same 
in  the  initial  and  final  states,  and  that  the  force  causing  the 
difference  in  momentum  flux  between  initial  and  final  states  is  the 

difference  in  pressures.  The  basis  of  these  equations,  and  their 

1  2 

limitations,  are  considered  in  standard  textbooks.  * 

It  is  convenient  to  work  with  these  equations  in  forms  which 
express  (1)  the  change  in  specific  enthalpy,  (h  -  hQ);  (2)  the  shock 
velocity,  ug,  into  stationary  initial  gas;  and  (3)  the  particle 
velocity,  u^,  imparted  by  passage  of  the  wave,  explicitly  in  terms 
of  the  initial  (subscript  zero)  and  final  (unsubscripted)  pressures, 
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P,  and  specific  volumes,  v.  Thus  we  have: 

(h  -  h  )  =  1/2(P  -  P  ) ( v  +  v) 
o  o  o 


u  =  [(P  -  P  )(v  -  v)]1/2 

p  o  o 


(1) 

(2) 

(3) 


In  addition  to  these  basic  equations,  we  use  the  thermal  and 
caloric  equations  of  state  which  relate  volume  and  enthalpy  to  tem¬ 
perature,  T,  and  pressure. 


t(v,P,T)  =  0 

(4) 

c(h,P,T)  =  0 

(5) 

For  ideal  gases  of  given  composition,  the  enthalpy  depends  only  on 
the  temperature.  For  this  reason,  temperature  is  used  as  the  inde¬ 
pendent  parameter  in  seeking  the  solution  of  equation  (1).  Shock 
wave  conditions  for  a  particular  value  of  some  parameter  such  as 
shock  velocity  or  particle  velocity  are  obtained  by  imposing  the 
desired  value  as  an  auxiliary  condition  and  varying  the  temperature 
for  which  equations  ( 1 ) — ( 5 )  are  solved  until  the  condition  is 
satisfied. 
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1.2  REFLECTED  SHOCK  WAVES 


Reflected  shock  waves  are  treated  by  equating  the  final  data 

from  the  original  shock,  P,  v,  and  h,  to  the  initial  conditions  for 

the  reflected  shock,  P^,  v^,  and  h^,  and  proceeding  to  solve  equations 

(l)-(5)  for  the  reflected  shock  parameters,  P',  v*,  h’,  u',  and  u'  for 

s  p 

an  assumed  T'.  The  coordinate  system  for  reckoning  velocities  in  the 
reflected  shock  wave  becomes  transformed  to  one  which  is  reversed  in 
direction  from  that  of  the  incident  shock  and  which  is  stationary,  not 
in  the  laboratory,  but  in  the  final  gas  behind  the  incident  shock. 

The  auxiliary  condition  which  matches  the  reflected  shock  to  the  in¬ 
cident  one,  and  thereby  fixes  T',  is  that  u*,  the  particle  velocity 
behind  the  reflected  shock  in  the  laboratory  system,  must  be  zero. 
Hence  this  condition  is: 


u*  =  u»  -  u  =o  (6) 

P  P  P 

The  reflected  shock  wave  velocity  in  the  laboratory  system  then 
becomes: 


u*  =  u'  -  u 
s  s  p 


1.3  CHEMICAL  REACTIONS 


(7) 


In  considering  shock  waves,  we  wish  to  restrict  our  attention  to 
solutions  of  equations  (l)-(5)  which  represent  compression  waves 
[(P  -  Pq)  >  0,  (vq  -  v)  >0]  with  velocities  greater  than  the  sound 
velocity.  Normally  we  are  also  interested  in  final  temperatures 
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which  are  greater  than  the  initial  temperature.  In  the  absence  of 
chemical  reactions  or  other  internal  sources  or  sinks  of  energy, 
these  conditions  are  equivalent.  However,  chemical  reactions  have 
profound  effects  on  the  equations  of  state  and  further  affect  the 
regimes  of  temperature  and  shock  velocity  in  which  solutions  of 
equation  (1)  represent  shock  waves. 

The  occurrence  of  endothermic  reaction  generally  raises  the 
shock  velocity  and  the  final  density  required  to  achieve  a  parti¬ 
cular  temperature,  so  that  no  complication  is  encountered  when  T  is 
kept  greater  than  T  .  It  can  happen,  though,  that  there  is  a  minimum 
shock  velocity,  greater  than  the  sound  velocity,  compatible  with  a 
given  extent  of  endothermic  reaction.  This  minimum  represents  ter¬ 
mination  of  solutions  of  equation  (1)  at  a  large  pressure  and  zero 
volume . 

When  exothermic  reactions  occur,  there  is  an  additional  compli¬ 
cation.  There  is  a  temperature,  T  9  greater  than  Tq,  which  represents 
reaction  at  constant  volume,  v  =  v^,  and  below  which  the  solution  of 
equation  (1)  requires  v  >  v^.  At  temperatures  above  T^,  the  solutions 
are  such  that  there  is  a  minimum  value  of  ug,  greater  than  the  sound 
velocity,  above  which  there  are  normally  two  P-v-T  conditions  for  a 
given  Ug  and  below  which  there  are  none. 

Of  all  the  stoichiometrically  possible  states  attainable  by 
chemical  reactions,  one  that  is  of  particular  interest  is  that  which 
satisfies  the  condition  of  chemical  equilibrium.  The  solutions  of 


equation  (1)  for  conditions  in  which  h  and  the  P-v-T  relationship 
are  governed  by  chemical  equilibrium  have  a  particular  locus  in  the 
P  -  v  plane,  the  so-called  equilibrium  Hugoniot  curve.  Shock  waves 
into  media  initially  at  equilibrium  are  simple  in  that  this  locus  in¬ 
cludes  the  initial  point,  and  weak  shocks  approach  no-reaction  behav¬ 
ior.  Shocks  into  metastable  initial  media,  with  the  final  state 
determined  by  chemical  equilibrium,  are  complicated  by  the  effects 
described  above  arising  from  exothermic  or  endothermic  reaction. 

1.4  CHAPMAN-JOUGUET  DETONATIONS 

The  Chapman- Jouguet  theory  treats  detonations  as  steady,  one¬ 
dimensional  shock  waves  in  which  the  final  state  is  reached  by 
exothermic  reaction.  .  The  detonation  properties  are  given  by  the 
solution  of  equations  ( 1 ) — ( 5 )  having  the  minimum  shock  velocity 

3 

compatible  with  equilibrium  in  the  chemically  reacting  system. 

This  solution  is  called  the  Chapman-Jouguet  solution,  and  is  deter¬ 
mined  by  the  tangency  condition: 

(dP/dv)  +  (P  -  P  )/(v  -  v)  =  0  (8) 

eq  Hug  o  o 

where  (dP/dv)^  ^Ug  is  the  slope  of  the  locus  of  equilibrium  solu- 

tions  of  equation  (1),  and  -  (P  -  P  )/(v  -  v)  is  the  slope  of  the 

o  o 

line  connecting  the  Hugoniot  point  P,v  with  the  initial  point,  ^0*v0» 
and  determines  ug  by  equation  (2). 


1.5  CHEMICAL  SYSTEMS 


r 


The  extensive  properties  of  the  medium  considered  in  equations 
(l)-(8)  are  on  a  specific  (per  gram)  basis.  The  properties  of  indi¬ 
vidual  chemical  species,  however,  are  most  conveniently  formulated  on 
a  molar  basis.  The  properties  of  the  system  may  then  be  formulated 
from  the  description  of  the  mixture  of  chemical  species. 

The  first  item  in  this  description  is  an  ordered  list  of  the 
chemical  species  included.  The  subscript  t_  is  used  to  denote  the 
position  of  each  species  in  this  list,  with  the  integer  s^  denoting 
the  length  of  the  list.  The  identity  of  the  t species  is  indicated 
by  its  chemical  formula.  The  composition  of  the  system  is  specified 
by  a  set  of  mole  fractions,  x^,  which  is  likewise  an  ordered  list 
(vector)  of  length  s_. 

Two  additional  sets  of  data,  likewise  ordered,  suffice  to  specify 
the  conditions  of  an  ideal  gas  mixture  at  a  given  temperature  and 
pressure.  These  are  the  molecular  weights,  M  and  the  molar  en¬ 
thalpies,  H^.  Thus  the  conditions  of  the  initial  gas  before  a  shock 
wave  are  given  in  terms  of  the  initial  (subscript  zero)  values  of 
these  parameters  by 


RT 

vo =  ftt 

o  o 


RT 


£ 

I 

t=i 


xt  Mt 
o 


(9) 


* 
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S  £ 

7  x  H.  V  x*  H 

.  t  t  .L,  t  t 

t=l  o  o  t=l  o  o 


y  x  m4 
+L  t  i 
t=l  o 


The  condition  of  the  system  after  chemical  reaction  is  described 
as  follows.  Let  nt  represent  the  numbers  of  moles  of  each  species 
formed  from  one  mole  of  the  original  gas.  Now  we  allow  the  possibil¬ 
ity  that  the  s**1  species  may  be  a  pure  condensed  (solid)  phase.  The 
control  integer  £  has  the  value  zero  if  the  s**1  species  is  gaseous, 
one  if  it  is  solid.  (If£=l,  xg  =0.)  We  then  formulate  the  mole 


number  sums 


S-£ 

J,  "t 


'solid  £-ns 


and  the  mole  fractions,  relative  to  the  gaseous  phase. 


t  n 


In  these  terms,  the  equations  of  state  become; 


n  RT  pn  V 
gas  A  s  s 

V  =  .  ...  —  --  4-  ■  i ,  -  

PM  M 

o  o 
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t 


(15) 


t=l 


XtHt 


where  V  is  the  molar  volume  of  the  solid  s**1  species. 

s  — 

1.6  CHEMICAL  EQUILIBRIUM 

4 

Chemical  equilibria  are  treated  by  the  Brinkley  method  as  ex¬ 
tended  to  two-phase  systems  by  Fickett.'*  No  outstanding  advantage  of 
the  Brinkley  method  over  other  available  methods  is  claimed;  the  choice 
of  this  method  has  been  based  on  local  familiarity  with  it.  The 
parameters  by  which  the  user  communicates  with  the  formalism  are  dis¬ 
cussed  here.  Additional  details  of  the  method  and  specific  examples 
of  its  use  are  provided  in  Appendix  A. 

The  formulas  of  the  chemical  species  are  represented  by  a 
matrix  whose  elements,  a  are  the  formula  coefficients  specifying 
the  number  of  atoms  of  each  element  in  each  species.  This  matrix  has 
a  row  (index  t)  for  each  of  the  s_  species  and  a  column  (index  k_)  for 
each  of  the  c'  chemical  elements. 

The  elemental  composition  of  the  system  is  specified  by  an 
ordered  set,  Q,  ,  of  length  c*,  whose  elements  are  the  numbers  of  gram 

K 

atoms  of  each  element  in  the  system.  From  the  previous  section,  the 

system  to  be  considered  in  a  shock  wave  calculation  is  one  mole  of 

initial  gas,  whose  composition  is  specified  by  .  Thus: 

o 


¥ 


v 


* 
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(16) 


* 


\ =  Jx  \°tk 


In  the  Brinkley  method t  chemical  equilibrium  and  stoichiometry 
are  actually  formulated  in  terms  of  a  matrix,  v£j*  anc^  a  exposition 
vector,  qy  derived  from  and  Q^.  This  derivation  requires  specifi¬ 
cation  of  a  set  of  c  stoichiometrically  independent  component  species 
(species  whose  formula  vectors,  a.  ,  form  a  c  by  c  matrix  whose 
determinant  is  not  zero).  This  number,  c_,  is  normally  equal  to  £f , 
the  number  of  chemical  elements.  It  may  be  smaller  than  c^  when  two 
or  more  elements  appear  only  in  fixed  combinations. 

This  specification  is  accomplished  by  an  auxiliary  vector  which 
assigns  a  unique  re-ordering  integer,  =  1,  . . .  to  each  species. 
The  selected  set  of  independent  species  are  assigned  elements  a^  =  1, 
...  £;  the  remaining  species,  for  which  a  =  (£  +  1),  ...  s,  are  then 
dependent.  Except  in  cases  where  the  s^  species  represents  a  second 
phase,  there  is  no  significance  to  the  ordering  within  these  two 
subsets. 

Free  energy  data  for  the  computation  of  equilibrium  constants 
are  communicated  as  a  set,  F^,  whose  elements  are  the  dimensionless 
standardized  chemical  potentials  for  each  species.  The  reference 
temperature  in  the  nondimensionalizing  factor,  RT,  is  the  temperature 
at  which  the  equilibrium  calculation  is  to  be  made.  The  basis  of  the 
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standardization  is  the  pressure  of  the  system,  P,  and  a  suitable  set 
of  standard  states  for  reckoning  enthalpies  of  formation. 


1.7  THERMODYNAMIC  DATA 

The  basis  for  reckoning  enthalpies  and  free  energies  for  all 
species  is  the  elements  in  their  standard  states  at  0°K.  Thus,  the 
molar  enthalpy  of  each  species  is  given  as 


Ht  =  <H  -  Vt  *  <AHf  >t 

O 


(17) 


where  (H  -  HQ)t  is  the  molar  enthalpy  of  the  t**1  species  relative  to 

itself  at  0°K,  and  (AHf  >t  is  the  molar  enthalpy  of  formation  of  the 

th  ° 

tx  species  from  the  elements  at  0°K. 

Because  the  range  of  initial  temperatures  normally  considered 
for  shock  waves  is  much  smaller  than  the  ranges  of  final  temperatures, 
and  because  free  energy  data  are  needed  only  for  the  final  temperature 
ranges,  different  functional  forms  are  used  to  express  (H  -  HQ)t  for 
the  different  ranges.  For  ideal  gases  over  the  range  of  initial  tem¬ 
peratures,  Tq,  the  following  polynomial  is  used: 


(H  -  Ho)t  =  at  +  btTo  +  ctT^  +  d^ 
o 


(18) 


For  final  temperature  ranges,  the  formulation  used  for  ideal  gases  is 


(H  -  Ho)t 
RT 


=  At  + 


BtT 


2 

C  T  + 


4 

DtT  +  EtT 


(19) 


* 
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The  corresponding  expression  for  the  standard  Gibbs  free  energy 
function  is  then: 


* 


% 


% 


% 


(F 


Vt 


RT 


=  At(l  -  In  T)  -  B  T  -  1/2C  T" 


i/3Dtr 


-  1/4E  T  -  k  (20) 


For  ideal  gases,  the  function  F*  referred  to  in  the  preceding 
section  takes  the  form: 


(F  -  H  )  <AHf  >t 

F*  =  - 2_L  +  - 2 —  +  in  P 

t  RT  RT 


(21) 


where  P  is  the  pressure  of  the  system,  in  units  of  the  standard 
pressure  (1  atmosphere)  for  reckoning  F°. 

When  a  solid  is  included  among  the  constituents  of  the  final 

state,  it  is  considered  to  have  a  constant  molar  volume,  V  ,  over  the 

’s’ 

limited  temperature  and  pressure  range  involved.  In  addition,  the 
small  effects  of  pressure  on  its  enthalpy  and  chemical  potential  are 
neglected.  Thus,  for  the  species,  equations  (19)  and  (20)  are 
used,  and  F^  is  given  by: 

S 

,  (F°  -  Vs  (tHf0>s 

Fs  =  RT - -  *  ~W~  *  <1  -  £J1"  P  (22> 


where  £  is  the  phase  control  integer  introduced  previously. 
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Chapter  2 


THE  HUG  CODE 


The  HUG  code  was  written  to  process  data  according  to  the 
equations  of  Chapter  1  using  the  IBM  7090  or  7094  computer  under  the 
MONITOR  system.  Two  of  its  major  sub-codes  are  written  in  FAP.  One 
of  these  is  the  code  for  solving  chemical  equilibria  by  the  method  of 

* 

Brinkley  as  extended  to  two-phase  systems  and  coded  by  Fickett. 

The  other  major  FAP  routine  is  the  linear  feedback  equation 

solver,  FROOT,  also  due  to  Fickett.6  This  routine  is  used  in  all  the 

shock  wave  computations  to  determine  the  pressure  which  satisfies 

equation  (1)  for  the  given  temperature.  In  addition,  it  is  used  to 

control  outer  iterations  to  determine  the  temperature  for  which  the 

auxiliary  conditions,  equations  (6),  (8),  or  an  arbitrarily  specified 

velocity,  u  or  u  ,  are  satisfied.  Each  use  of  the  FROOT  routine 

s  p 

requires  the  specification  of  a  value,  e,  such  that  a  value  of  the 
variable  satisfying  the  equation  with  a  residual,  6  ,  is  accepted  as 
a  solution  when  6<e  . 
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The  remainder  of  the  HUG  code  is  written  in  FORTRAN  II.  A  great 
deal  of  the  code  is  logical  control  which  governs  the  many  options 
available  to  the  user.  A  run,  or  job,  of  computing  may  consist  of 
one  or  many  separate  problems,  each  of  which  may  contain  one  or 
several  individual  calculations,  or  points.  Most  of  the  options  are 
selected  by  means  of  a  set  of  control  integers,  which  are  specified 
anew  for  each  problem  within  the  run.  Data  on  the  chemical  system 
may  be  respecified  (block  by  block)  between  problems,  but  if  no  change 
is  to  be  made  from  the  previous  problem,  most  of  the  data  are  re¬ 
tained  and  re-used  as  called  for  in  subsequent  problems.  Likewise, 
data  specifying  the  conditions  for  the  sequence  of  points  in  the 
problem  may  be  changed  between  problems . 

Provision  is  also  made  for  the  optional  loading  of  certain  prob¬ 
lem  data  which  would  otherwise  be  computed  internally.  These  data, 
if  loaded,  are  used  in  preference  to  the  results  of  the  internal 
computations.  They  are  destroyed  at  the  end  of  the  problem  for  which 
they  were  loaded,  and  hence  must  be  loaded  anew  for  each  problem  in 
which  they  are  to  be  used. 

2.1  INPUT  SPECIFICATIONS 

Input  data  are  entered  by  the  instruction:  READ  INPUT  TAPE  IN, 
etc.  The  value  of  IN  is  assigned  in  the  main  program.  At  LASL,  the 


value  10  is  used. 


The  various  types  of  input  data  used  by  the  code,  and  the  for¬ 
mats  for  the  cards  from  which  they  are  transcribed  onto  the  input 
tape,  are  specified  below.  All  data  prescribed  below  are  decimal 
numbers.  Table  I,  at  the  end  of  this  section,  summarizes  the  loading 
formats. 

I .  LABEL  CARD 

This  card  carries  labels  for  certain  columns  in  the  printout  and 
values  of  e,  the  acceptable  limit  of  proximity  to  zero,  for  the  four 
FROOT  iterations.  It  is  the  first  data  card  loaded  in  a  run,  and  it 
is  loaded  only  once. 

Columns  1-30,  Labels  (_  Formats  A3) 

Hollerith  characters  are  entered  in  ten  fields  of  width  three. 
Those  in  the  first  six  fields  label  the  columns  (chemical  elements)  in 
the  a-matrix.  (See  a-MATRIX  CARDS,  columns  61-72).  The  label  con¬ 
vention:  C,  H,  0,  N,  M  (noble  gas),  X  (halogen)  is  used  at  LASL.  The 
contents  of  the  last  four  fields  label  the  a-vectors. 

Columns  31-70,  Definitions  of  Zero  (Formats  E10.2) 

Positive  floating  point  numbers  containing  the  decimal  point  with 
one  digit  to  the  left  and  two  digits  to  the  right  of  it,  (+X.XX*EE), 
are  entered  in  four  fields  of  width  ten. 

Columns  33-40,  ep 

In  the  iteration  on  pressure  to  solve  equation  (1)  for  a 
specified  temperature,  a  value  satisfying 
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♦ 


(P  -  P  )  ( V  +  v) 
o  o 

82. 58 56 ( h  -  h  ) 
o 


is  acceptable  as  a  solution.  P  is  pressure  (atmospheres),  v  is 
specific  volume  (cc/gram),  and  h  is  specific  enthalpy  (calories/gram). 
A  value  of  of  +1.0  x  10  ^  usually  gives  sufficient  accuracy. 

Columns  43-50, 

In  the  iteration  on  temperature  to  achieve  a  specified  shock 
wave  or  particle  velocity,  the  solution  is  accepted  when 


1.0  - 


DSTORE 


=  6  J  <  e 


T* 


where  u  is  the  velocity,  ug  or  u^, 
prespecified  value.  A  value  of 
machine  accuracy. 

Columns  53-60,  e 

K 


being  tested  and  DSTORE  is  the 
of  +1.0  x  10  6  provides  near 


The  solution  for  an  assumed  temperature  for  a  shock  wave  prop¬ 
agating  into  the  medium  behind  a  given  incident  shock  wave  is 
accepted  as  the  reflected  shock  solution  when 

I  u*  i 


1.0  -  -2- 
u 

p 


<  e. 


9 


where  u^  is  the  particle  velocity  behind  the  incident  shock  (ambient 
gas  stationary)  and  uf  is  the  particle  velocity  (coordinate  system 
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* 


fixed  in  the  gas  behind  the  incident  shock)  being  tested.  A  value  of 

“6 

e  of  +1.0  x  10  provides  near  machine  accuracy. 

R 

Columns  63-70, 

The  shock  wave  solution  for  an  assumed  temperature,  T,  is 
accepted  as  a  Chapman-Jouguet  detonation  solution  when 

(v  -  v) 

p"7  ( AP/Av)  +  1.0 
o 

where  P  and  v  are  pressure  and  specific  volume  at  the  temperature  T, 

and  (AP/Av)  is  the  approximate  slope  of  the  reactive  Hugoniot  curve 

(equilibrium,  or,  if  desired,  fixed  exothermic  composition)  obtained 

by  evaluating  AP  and  Av  as  the  differences  between  P  and  v  at  two 

Hugoniot  points  displaced  from  T  by  +  6T  and  -  6T •  The  size  of  6T 

is  chosen  internally  as  a  minimum  needed  to  make  AP/P  >  0.001  and  Av/v 

>  0.001,  thereby  assuring  numerically  meaningful  values.  A  value 
-4 

of  Cp  of  +1.0  x  10  usually  provides  near  machine  accuracy. 

Lack  of  precision  may  prevent  acceptable  convergence  of  these 
iterations  under  unfavorable  circumstances.  When  cannot  be  satis¬ 
fied,  owing  to  near  indeterminacy  when  (h  -  hQ)  ~  0  or  other  reasons, 
the  working  value  of  T  is  increased.  Precision  in  this  solution  must 
be  maintained  to  avoid  compounded  difficulties  in  the  outer  iter¬ 
ations.  Provision  is  made  for  temporary  increases  in  e^,  eR,  and  eD 
when  convergence  cannot  be  obtained  with  the  original  values* 
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Fivefold  increases  are  made,  and  notes  are  printed  out  when  this  is 
done.  If  convergence  still  fails  after  two  increases,  the  velocity 
or  reflected  shock  point,  or  the  entire  detonation  problem,  is 
rejected.  The  original  e*s  are  restored  when  this  occurs.  They  are 
also  restored  before  all  new  problems. 

II.  MPROB  CARD  (Formats  13) 

This  card  carries  the  problem  control  integers  which  govern  the 
principal  options  in  loading,  execution,  and  printing  in  each  problem. 
It  must  be  loaded  following  the  LABEL  card  in  the  first  problem  of  a 
run,  and  as  the  first  data  card  in  subsequent  problems.  All  entries  on 
this  card  are  unsigned,  fixed-point  integers,  right  shifted  in  fields 
of  width  three.  Eight  such  fields,  beginning  with  column  1,  and  five 
additional  fields,  beginning  with  column  37,  are  used.  The  following 
list  gives  the  numbers  to  be  supplied,  their  possible  values,  and  the 
significance  of  each  value. 

Column  3 

MPRO  Specifies  the  type  of  calculation  to  be  made. 

1  Shock  wave(s)  in  the  initial  medium. 

2  Detonation  wave(s)  in  the  initial  medium. 

3  Equilibrium  calculation( s)  only,  no  shock  waves. 

Column  6 

MINC  Specifies  the  chemical  change  on  the  incident  shock 
wave. 

1  No  chemical  reaction. 
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2  Final  composition  specified  by  data  loaded  on  type  03 
general  data  card(s). 

3  Equilibrium  final  composition. 

Column  9 

MREF  when  MPRO  =  1  or  2  specifies  the  occurrence  of,  and 

chemical  change  on,  reflected  shock  waves  following 
each  incident  shock  calculation. 

0  No  reflected  shock  calculations. 

1  Reflected  shock  wave  with  composition  the  same  as  that 
behind  the  incident  shock. 

2  Reflected  shock  with  composition  specified  by  data 
loaded  on  type  03  general  data  card(s). 

3  Reflected  shock  with  equilibrium  final  composition. 

4  Reflected  shock  with  equilibrium  final  composition; 

independent  components  for  equilibrium  calculation 
specified  differently  for  incident  and  reflected 
shocks.  (Usable  only  if  MINC  =  3,  0  =  0) 

MREF  when  MPRO  =  3  specifies  a  set  of  pressures  for  equilibrium 
calculations.  A  set,  even  if  it  be  of  length  one,  must 
be  specified. 

1  Arithmetic  progression  of  pressures,  specified  by  a 
type  08  general  data  card. 

2  Sequence  of  pressures,  P,  individually  specified  by 
type  09  general  data  card(s). 

Column  12 

MAVECT  Specifies  the  number  of  a-vectors  available  in  storage, 
0,  1,  2,  3,  or  4.  This  number  is  used  whenever 
equilibrium  calculations  are  called  for. 
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Column  15 


MITER  when  MPRO  =  1,  specifies  whether  shock  wave  calculations 
are  to  be  made  for  a  sequence  of  final  temperatures, 
shock  velocities,  or  particle  velocities.  A  sequence, 
even  if  of  length  one,  must  be  specified. 

Care  must  be  used  in  specifying  such  series  for  exo¬ 
thermically  reactive  mixtures.  Temperatures  for  which 
v>vq  are  automatically  increased  five  percent,  as  many 
as  fifty  times  if  need  be,  until  a  solution  is  obtained 

for  u  .  Shock  velocities  for  which  a  solution  cannot  be 
s 

found  are  rejected,  and  the  next  value  in  the  sequence 

is  tried.  When  two  solutions  exist  for  a  given  u  , 

s 

either  one  may  be  obtained,  although  the  higher  tem¬ 
perature  one  is  favored  by  the  code’s  trial  values. 

1  Arithmetic  progression  of  final  temperatures,  T, 
specified  by  a  type  06  general  data  card. 

2  Sequence  of  final  temperatures,  T,  individually 
specified  by  type  07  general  data  card(s). 

3  Arithmetic  progression  of  shock  velocities,  u  , 
specified  by  a  type  08  general  data  card. 

4  Sequence  of  shock  velocities,  u  ,  individually 
specified  by  type  09  general  dafa  card(s). 

5  Arithmetic  progression  of  particle  velocities,  u  , 

specified  by  a  type  08  general  data  card.  ^ 

6  Sequence  of  particle  velocities,  u  ,  individually 
specified  by  type  09  general  data  Bard(s). 
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MITER  when  MPRO  =  2,  specifies  the  choice  of  a  set  of  initial 
conditions,  temperatures  or  pressures,  for  a  sequence 
of  detonation  calculations.  A  sequence,  even  if  of 
length  one,  must  be  specified. 

1  Arithmetic  progression  of  initial  temperatures,  Tq, 
specified  by  a  type  06  general  data  card, 

2  Sequence  of  initial  temperatures,  Tq,  individually 
specified  by  type  07  general  data  card(s). 

3  Arithmetic  progression  of  initial  pressures,  PQ, 
specified  by  a  type  08  general  data  card. 

4  Sequence  of  initial  pressures,  Pq,  individually 
specified  by  type  09  general  data  card(s). 

MITER  when  MPRO  =  3,  specifies  the  choice  of  a  set  of  tem¬ 
peratures,  T,  for  a  sequence  of  equilibrium  calcu¬ 
lations.  A  sequence,  even  if  of  length  one,  must  be 
specified. 

1  Arithmetic  progression  of  temperatures,  T,  specified 
by  a  type  06  general  data  card. 

2  Sequence  of  temperatures,  T,  individually  specified  by 
type  07  general  data  card(s). 

Column  18 

MCONC  controls  optional  printing  of  molar  concentrations,  £ 
n 

C  =  x  for  shock  and  detonation  problems. 

— t  t  M  v 

o 

(MPRO  =  1  or  2.)  Solid,  if  present,  is  treated  as 
dispersed  in  the  gas. 

P/RT  for  equilibrium  problems  (MPRO  =  3). 
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In  this  case,  concentrations  are  based  on  the  gas 
volume,  and  if  £  =  1,  should  be  ignored. 

0  Concentrations  not  printed. 

1  Concentrations  printed. 

Column  21 

MALPHA  Controls  loading  of  a-matrix  cards. 

0  No  a-matrix  cards  to  be  loaded. 

1  £  a-matrix  cards  to  be  loaded. 

Column  24 

MA  Specifies  the  number  of  a-vectors,  0,  1,  2,  3,  or  4,  to 

be  loaded  on  ^-vector  card(s). 

Column  39 

c_  The  number  of  independent  components  for  equilibrium 

computations,  which  must  also  equal  the  number  of 
nonzero  columns  in  the  a-matrix.  (See  Appendix  A, 
Example  3.)  £<6. 

Columns  40-42 

£  The  number  of  species  in  the  system  for  shock  wave 

and/or  equilibrium  computations.  £<100 

Column  45 

£  Phase  control  integer  discussed  in  Section  1.5.  The 

value  of  £  loaded  here  is  used  when  0=0. 

0  One  phase  present,  a  mixture  of  ideal  gases. 


1  Two  phases  present,  second  phase  a  pure  species.  A 

negative  value  of  x  may  be  obtained,  indicating  an 
artificial  solution. 

Column  48 

p*  Trial  value  for  £,  used  for  first  calculation  when 

0=1.  Subsequent  calculations  use  the  value  of  £  from 
the  previous  calculation  as  the  trial  value, 

0  One-phase  system  tried  first. 

1  Two-phase  system  tried  first. 

Column  51 

0  Controls  use  of  the  evaporation-saturation  test  to 

determine  £. 

0  £  used  as  loaded. 

1  £  determined  by  code. 

When  this  option  is  used,  the  code  must  have  at  least 
one  equilibrium  equation  to  consider  in  addition  to  the 
question  of  saturation;  thus  £  must  be  at  least  three. 

If  £  =  0  is  determined,  the  printed  values  of  xg  and  C^ 
are  artificial  and  should  be  negligibly  small. 

III.  a-MATRIX  CARDS 

These  cards  carry  the  chemical  symbol,  the  molecular  weight, 
enthalpy  data  for  a  range  of  initial  temperatures,  and  the  chemical 
formula  of  each  species.  A  system  of  £  species  requires  £  such  cards. 
These  must  follow  the  MPROB  card  and  be  loaded  when  MALPHA  =  1.  The 
order  of  species  used  for  loading  and  printing  data  throughout  the 
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problem  is  determined  by  the  order  of  the  n-matrix  cards.  Data 
loaded  by  these  cards  are  retained  and  are  available  for  use  in  sub¬ 
sequent  problems.  The  entries  on  the  t**1  card  are  given  below: 
Columns  1-6,  Species  Symbol  (Format  A6) 

Hollerith  characters  loaded  here  are  printed  as  the  tth  entry 
in  the  ordered  list  of  species  in  the  system,  and  serve  to  label  the 
rows  of  other  data  for  that  species. 

Columns  7-12,  Molecular  Weight  (Format  F6.3) 

The  molecular  weight  of  the  t**1  species  is  entered  as  a  fixed- 
point  number  with  the  decimal  point  assumed  to  be  between  the  third 
and  fourth  positions  in  the  field.  Zeros  to  the  right  of  the  decimal 
point  should  be  entered.  The  molecular  weight  is  used  in  shock  and 
detonation  wave  problems  to  compute  Mq,  as  defined  in  equation  (9). 

It  is  used  only  for  those  species  with  nonzero  x^  ,  and  need  not  be 

o 

entered  for  other  species,  if  the  user  so  chooses. 

Columns  13-60,  Initial  Enthalpy  Coefficients  (Formats  E12.7) 

The  polynomial  coefficients  for  the  molar  enthalpy  of  the  t 
species  in  the  range  of  initial  temperatures  for  shock  and  detonation 
wave  problems  are  entered  in  four  fields  of  width  twelve.  These 
coefficients  are  defined  in  equation  (18),  and  are  entered  in  the 
order  d^,  c^,  bt,  at,  as  signed,  eight-digit,  floating-point  numbers 
with  the  decimal  point  assumed  to  be  to  the  right  of  the  first  digit. 
The  unit  of  enthalpy  in  equation  (18)  is  calories  per  gram  mole;  the 
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unit  of  temperature  is  °K.  For  species  which  will  not  be  included  in 
the  initial  gas  mixture,  these  fields  may  be  left  blank. 

Columns  61-72,  The  q-matrix  (Formats  12) 

Six  fields  of  width  two  are  provided  for  entering  the  coefficients 
of  up  to  c  5  6  independent  elements  in  the  chemical  formula  for  the 
t**1  species,  i.e.,  the  t^  row  of  the  a-matrix.  These  data  are  needed 
for  all  equilibrium  computations.  The  element  represented  by  each 
field  is  fixed  by  the  user,  and  should  correspond  to  the  labelling  in 
columns  1-18  of  the  LABEL  card.  When  fewer  than  six  elements  are 
involved  in  the  system,  any  c_  fields,  not  necessarily  consecutive, 
may  be  selected  for  use.  The  remaining  (6  -  c.)  fields  may  be  zero  or 
blank  for  all  species;  such  superfluous  columns  are  automatically 
omitted  from  the  printed  output.  The  numbers  in  the  a-matrix  are 
entered  as  fixed-point  integers,  right  shifted  in  their  fields. 

Zeros  should  be  punched  to  avoid  -0's  appearing  in  the  printout. 
Positive  signs  need  not  be  punched.  Negative  signs  may  be  used;  they 
are  useful  in  ionization  and  restricted  equilibrium  problems,  as  dis¬ 
cussed  in  Appendix  A. 

IV.  a-VECTOR  CARDS  (Formats  13) 

These  cards  carry  the  elements,  a.t,  of  the  renumbering  vectors 
by  which  are  designated  the  independent  component  species  for  equilib¬ 
rium  computations.  (See  Section  1.6.)  Each  a-vector  consists  of  s_ 
ordered  positive  integers,  which  are  loaded  as  unsigned  fixed-point 
numbers,  right  shifted  in  up  to  twenty-four  consecutive  fields  of 
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width  three  per  card.  Each  £-vector  begins  with  the  first  field  of  a 
card  and  when  £  <  24,  is  completely  contained  on  that  card.  Unused 
fields  at  the  right  of  a  card  are  left  blank.  Loading  is  controlled 
by  the  values  of  £  and  MA,  and  the  block  of  a-vectors  must  be  loaded 
following  any  a-MATRIX  cards  whenever  HA  >  0.  Data  so  loaded  are 
retained  for  subsequent  problems  if  not  replaced. 

For  a  given  set  of  £  chemical  species  composed  of  £  independent 
elements,  any  of  several  selections  of  component  species  which  will 
satisfy  the  requirement  of  independence  are  normally  possible.  As  a 
practical  requirement  of  the  code,  however,  these  must  be  selected  so 
that  their  mole  fractions  in  the  equilibrium  solution  will  exceed  a 
minimum  value,  presently  set  at  1.0  x  10~^.  This  requirement  often 
restricts  the  usefulness  of  £-vectors  to  particular  ranges  of  tem¬ 
perature,  pressure,  or  composition  conditions,  and  occasionally  pre¬ 
cludes  equilibrium  calculations  in  precisely  stoichiometric  mixtures 
at  low  temperatures. 

For  this  reason,  provision  is  made  for  the  user  to  supply  alter¬ 
nate  a-vectors,  to  be  used  should  the  original  selection  prove  in¬ 
adequate  during  the  course  of  the  problem.  Space  is  provided  for  as 
many  as  four  a-vectors.  For  most  problems  (MREF  <4,  0  =  0)  only  one 
is  needed,  and  only  two  are  available.  For  MREF  =  4  or  for  0=1, 
two  £-vectors  are  needed,  and  alternates  for  each  are  available.  The 
selection  of  £-vectors  under  the  various  options  is  summarized  in 
tabular  form  on  the  following  page. 
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a- VECTOR  SELECTION 
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For  two-phase  systems,  there  are  the  following  additional  re¬ 
quirements  on  the  ^-vectors.  When  0  =  o,  £  =  1,  it  is  required  that 
all  a  =  c.  When  0  -  1,  it  is  required  that  al  =  c,  a2  =  s,  and,  if 

the  alternate  set  be  used,  a3  =  c,  a4  =  s. 

—  s  —  —  s 

V.  GENERAL  DATA  CARDS  (Formats  7E10.5,  12) 

All  data  not  prescribed  in  the  foregoing  special  formats  are 
loaded  in  a  general  floating-point  format.  These  general  data  cards 
are  loaded  following  the  special  format  cards,  and  a  variable  number 
of  such  cards  may  be  included  with  each  problem.  Data  are  loaded  in 
the  general  data  format  as  signed,  six-digit,  floating-point  numbers 
with  the  decimal  point  assumed  to  be  to  the  right  of  the  first  digit, 
viz.,  ‘XXXXXX^E.  Fields  for  seven  such  numbers  occupy  columns  1-70, 
and  the  type  of  data  being  loaded  is  prescribed  by  a  load  control 
integer,  an  unsigned  fixed-point  number  having  a  value  between  one  and 
ten,  punched  in  columns  71-72.  Data  are  loaded  field  by  field.  When 
a  blank  field  is  encountered,  no  further  data  are  read  from  that  card; 
when  a  card  with  the  first  field  blank  is  encountered,  no  further  data 
are  read  for  that  problem. 

The  data  associated  with  each  load  control  integer  are  enumerated 
below: 

TYPE  01,  THERMO 

These  cards,  which  contain  the  digits  01  in  columns  71-72,  carry 
the  constants  from  which  the  molar  thermodynamic  functions  for  each 
species  are  computed  by  equations  (17),  (19),  and  (20)  in  the 
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temperature  range  of  the  equilibrium  and  final  shock  solutions.  The 

unit  of  enthalpy  in  equation  (17)  is  calories  per  gram  mole.  The  left 

hand  sides  of  equations  (19)  and  (20)  are  dimensionless;  the  unit  of 

temperature  on  the  right  hand  sides  is  °K.  A  system  of  s_  species  uses 

data  from  s  type  01  cards;  the  order  of  cards  must  match  that  of  the 

o -MATRIX  cards.  The  t**1  such  card  contains,  in  order,  the  quantities 

Et,  Dt,  C  ,  Bt,  At,  kt,  and  (AHf  )t»  for  the  Jtth  species.  These  data 

o 

are  retained  for  use  in  subsequent  problems  if  not  replaced.  Since 
only  one  set  of  these  quantities  is  available  in  a  problem,  its  range 
of  validity  must  include  the  temperatures  of  both  incident  and  re¬ 
flected  shock  waves,  when  reflected  shocks  are  being  considered. 

TYPE  02,  XCOMPO 

On  these  cards,  which  contain  the  digits  02  in  columns  71-72,  are 

loaded  the  ordered  list  of  initial  mole  fractions,  xt  ,  for  use  in 

o 

shock  and  detonation  wave  problems.  For  equilibrium  problems, 

(MPRO  =3),  the  formal  composition  of  the  system  may  be  loaded  in 

this  way,  or  use  may  be  made  of  the  type  10  card.  The  sum  of  the  s. 

values  of  x  must  be  unity,  to  within  one  part  in  10^,  or  problems 
o 

with  MPRO  < 3  will  be  rejected.  This  block  of  data  is  retained  for 
subsequent  problems  if  not  replaced.  Clearly,  zeros  must  be  loaded 
when  they  occur  between  nonzero  values  in  the  list;  zeros  at  the  end 
of  the  list  need  not  be  loaded,  unless  they  should  be  needed  to  over¬ 
write  values  retained  from  a  previous  problem.  Thus  when  s_  is  large. 
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it  is  well  to  order  the  species  so  that  those  having  nonzero  x^ 

o 

appear  early, 

TYPE  03,  XMOLE 

These  cards,  which  contain  the  digits  03  in  columns  71-72,  are 

used  to  load  the  final  composition  in  shock  and  detonation  problems 

having  MINC  and/or  MREF  equal  to  2.  The  quantities  loaded  are  the 

ordered  list  of  mole  numbers,  n^,  defined  in  Section  1.5,  These  data 

are  retained  for  use  in  subsequent  problems,  and  the  same  precautions 

apply  as  in  the  case  of  XCOMPO  data. 

The  numbers  n^  are  readily  generated  for  arbitrary  extents  of 

specific  chemical  reactions  by  equating  the  original  mole  fractions, 

x^  ,  to  original  mole  numbers,  n^  ,  and  adding  or  subtracting  from 
o  o 

these  numbers  by  the  ordinary  rules  of  stoichiometry.  The  resulting 
values  of  n^  will  not,  in  general,  sum  to  unity,  but  they  will  re¬ 
present  the  proper  mass  of  material,  namely  Mq  grams.  Equation  (16) 
remains  satisfied  when  a  proper  set  of  n^  values  are  substituted  for 

Xt  ‘ 
o 

TYPE  04,  XCOMP 

These  cards,  which  contain  the  digits  04  in  columns  71-72,  are 
used  to  load  estimates  of  the  mole  fractions,  x^,  to  be  used  in  the 
initial  cycle  of  the  iteration  in  chemical  equilibrium  computations. 
Such  data  are  not  retained;  rather,  the  results  of  the  preceding 
computation  are  used  to  begin  the  iteration  in  each  succeeding 


37 


computation.  On  entry  to  each  new  problem,  all  x^  estimates  are  set 
to  0.1,  which  value  may  be  replaced,  in  sequence,  by  values  loaded 
by  type  04  general  data  cards.  In  fact,  only  the  estimated  mole 
fractions  for  the  independent  component  species  are  used,  the  remain¬ 
ing  values,  whether  0.1  or  some  value  specified  by  the  user,  are 
superfluous.  In  practice,  most  problems  are  solved  successfully, 
though  perhaps  less  efficiently,  with  the  initial  estimates  of  the 
x's  of  0.1,  and  use  of  the  type  04  card  is  seldom  necessary. 

TYPE  05,  CONDO 

This  card,  with  the  digits  05  in  columns  71-72,  is  used  to  load 
a  collection  of  data  specifying  initial  conditions,  the  molar  volume 
of  a  solid  sth  species,  and  a  trial  value  of  temperature  to  be  used 
in  detonation  and  certain  shock  wave  calculations. 

In  the  first  field  is  loaded  the  initial  temperature,  T  ,  in  °K, 
for  all  shock  problems  (MPRO  =  1)  and  those  detonation  problems 
(MPR0  =2)  with  MITER  >  2.  The  value  loaded  is  retained  for  sub¬ 
sequent  problems  if  not  replaced. 

In  the  second  field  is  loaded  the  initial  pressure,  P  ,  in 
atmospheres,  for  all  shock  problems  (MPRO  =  1)  and  those  detonation 
problems  (MPRO  =  2)  having  MITER  <  3.  The  value  loaded  is  retained 
for  subsequent  problems  if  not  replaced. 

In  the  third  field  may  be  loaded  the  value  of  the  specific 
volume,  vq,  in  cc/gram,  of  the  initial  medium  for  a  shock  wave  or 
detonation  problem.  The  value  of  vq  is  normally  computed  internally 
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by  equation  (9),  and  this  is  done  unless  a  value  other  than  zero  is 
loaded  by  the  user.  For  special  problems  involving  initial  media 
other  than  mixtures  of  ideal  gases,  vq  may  be  obtained  externally  and 
loaded.  A  value  so  loaded  is  used  only  in  the  problem  for  which  it 
is  loaded,  and  zeroed  before  a  new  problem  is  entered. 

In  the  fourth  field  may  be  loaded  the  value  of  the  specific 
enthalpy,  hQ,  in  calories  per  gram,  of  the  initial  medium  for  a  shock 
or  detonation  wave  problem.  The  value  of  hQ  is  normally  computed 
internally  by  equations  (10), (17),  and  (18),  and  this  is  done  unless 
a  value  other  than  zero  is  loaded  by  the  user.  This  loading  option 
is  useful  for  nonideal  gas  mixtures  and  for  shocks  into  media  outside 
the  normal  range  of  ambient  temperatures  where  the  coefficients  in 
equation  (18)  are  valid,  e.g.,  the  gas  behind  another  shock  wave.  A 
loaded  value  of  hQ  is  used  only  in  the  problem  for  which  it  is  loaded, 
and  is  zeroed  before  a  new  problem  is  entered. 

.  "t  h 

In  the  fifth  field  is  loaded  the  molar  volume  of  a  solid  s_ 

species,  V  ,  in  cc/gram  mole.  This  value  is  used  only  for  two-phase 
s 

problems;  it  is  retained  for  use  if  called  for  in  subsequent 
problems. 

In  the  sixth  field  may  be  loaded  an  estimated  value  of  the  final 

temperature,  T,  in  °K,  for  detonation  problems  (MPRO  =  2)  or  shock 

problems  (MPRO  =  1)  in  which  a  particular  velocity  condition  is 

specified  (MITER  >2).  If  no  estimate  is  loaded,  a  value  of 

T  =  10T  is  used  for  the  first  trial.  The  first  value  is  altered  by 
o  J 
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100°K,  in  the  appropriate  direction,  for  the  second  trial  before  the 
linear  feedback  routine  is  entered.  For  some  problems,  such  as  deto¬ 
nations  in  highly  diluted  mixtures,  these  values  may  be  inappropriate, 
and  the  user  may  need  to  supply  a  more  appropriate  trial  value.  A 
loaded  value  is  not  retained,  but  is  zeroed  before  a  new  problem  is 
entered. 

TYPE  06,  TEMPR 

This  card,  containing  the  digits  06  in  columns  71-72,  is  used  to 
load  data  specifying  an  arithmetic  progression  of  temperatures, 

Ti  =  T1  +  (i  -  1)  AT’  i  =  1*  ***  N6*  The  first  field  contains  Ng, 
an  integer  loaded  as  a  floating-point  number.  The  second  field  con¬ 
tains  T^,  in  °K;  the  third  field  contains  AT,  likewise  in  °K,  which 
may  be  positive  or  negative.  Values  loaded  by  this  card  are  retained 
for  subsequent  problems  if  not  replaced. 

It  is  often  well,  in  cases  where  shock  wave  solutions  may  not 
exist  at  low  temperatures,  or  where  the  equilibrium  computation  may 
encounter  difficulties  at  the  lowest  temperatures,  to  run  the  sequence 
in  order  of  decreasing  temperature. 

TYPE  07,  TEMPS 

These  cards,  which  contain  the  digits  07  in  columns  71-72,  are 
used  to  specify  an  arbitrary  sequence  of  as  many  as  twenty  temper¬ 
atures,  T^,  i  =  1,  ...  N?.  The  first  field  contains  the  value  of  N?, 
an  integer  loaded  as  a  floating-point  number.  Succeeding  fields 
contain  the  temperatures,  in  °K.  Values  loaded  on  these  cards  are 
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retained  for  subsequent  problems. 

TYPE  08,  VELR 

This  card,  which  contains  the  digits  08  in  columns  71-72,  is 
used  to  load  data  specifying  an  arithmetic  progression  of  velocities 
or  pressures, 

u.  *  u.  +  (i  -  l)Au,  or  P.  =  P.  +  (i  -  1)AP,  i  =  1,  ...  N  . 

IX  IX  o 

The  first  field  contains  the  number  N  ,  an  integer  loaded  as  a 

O 

floating-point  number.  The  second  field  contains  u^,  in  cm/sec,  or 
P^,  in  atmospheres.  The  third  field  contains  the  value  of  Au  or  AP 
in  the  above  units,  and  may  be  positive  or  negative.  Data  loaded  by 
this  card  are  retained  for  subsequent  problems. 

TYPE  09,  VELS 

These  cards,  which  contain  the  digits  09  in  columns  71-72,  are 
used  to  specify  an  arbitrary  sequence  of  as  many  as  twenty  velocities 
or  pressures,  u.  or  P^,  i  =  1,  ...  N^.  The  first  field  contains  the 
number  N^,  an  integer  loaded  as  a  floating-point  number.  The  succeed¬ 
ing  fields  contain  the  individual  values  of  u,  in  cm/sec,  or  P,  in 
atmospheres.  Data  loaded  by  these  cards  are  retained  for  subsequent 
problems  if  not  replaced. 

TYPE  10,  QATOM 

This  card,  which  contains  the  digits  10  in  columns  71-72,  may  be 
used  to  load  the  elemental  composition  of  the  system,  Q^,  i  =  1, 

...  c,  for  use  in  equilibrium  computations.  This  is  done  by  loading 
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in  c  <  6  consecutive  fields  the  values  of  Q  for  the  chemical  elements 

having  nonzero  columns  of  the  a-matrix.  Values  so  loaded  are  used  for 

all  equilibrium  calculations  within  the  problem  for  which  they  are 

loaded,  and  are  discarded  before  subsequent  problems  are  entered. 

Normally  the  are  computed  internally  from  the  initial  mole 

fractions,  ,  and  the  ex-matrix  by  equation  (16),  in  order  that  the 
o 

initial  and  final  volumes  and  enthalpies  computed  for  shock  problems 
be  for  the  same  mass  of  material.  However,  for  pure  equilibrium  prob¬ 
lems,  (MPRO  =  3),  the  initial  conditions  are  not  relevant,  and  for 
these  problems  and  for  certain  restricted  equilibrium  problems  dis¬ 
cussed  in  Appendix  A,  the  user  may  wish  to  load  the 
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TABLE  I  (continued) 

80  COLUMN  ENTRY 


44 


TABLE  I  (continued) 

80  COLUMN  ENTRY 
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2.2  OUTPUT  SPECIFICATIONS 

Data  to  be  printed  as  output  are  processed  by  the  instructions 
WRITE  OUTPUT  TAPE  10,  etc.  The  value  of  10  is  assigned  in  the  main 
program.  The  value  9  is  used  at  LASL. 

The  initial  page  of  output  for  each  problem  contains  the  input 
and  control  data,  appropriately  labelled,  and  the  results  of  certain 
preliminary  computations.  The  numbers  of  gram  atoms  of  the  elements, 
Q.,  are  printed  in  the  order  of  the  nonzero  columns  of  the  a -matrix. 
QBAR  is  a  vector  expressing  the  composition  of  the  system  as  formal 
mole  fractions  of  the  independent  component  species.  NUMAT  is  a 
matrix  similar  to  the  a -matrix  which  expresses  by  formula  coefficients 
the  chemical  formulas  of  the  dependent  components  in  terms  of  those 
of  the  independent  ones.  QBAR  and  NUMAT  are  computed  internally 
from  the  composition,  the  ct— matrix,  and  the  a— vectors. 

The  results  of  shock  or  detonation  wave  problems,  reflected  shock 
problems,  and  equilibrium  problems  are  printed  as  unlabelled  blocks 
of  data  for  each  point  in  the  problem.  The  formats  of  these  blocks 
are  given  on  the  following  pages • 
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SHOCK  AND  DETONATION  CALCULATIONS 
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REFLECTED  SHOCK  CALCULATIONS 
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2.3  OPERATION 


Assuming  proper  specification  of  input,  there  should  be  no 
machine  stops.  For  many  errors  in  input  and  for  some  other  diffi¬ 
culties,  the  code  will  write  a  note  briefly  describing  the  difficulty 
and  then  proceed  to  the  next  problem.  A  machine  stop  will  occur  if  a 
supposedly  independent  set  of  components  are  not,  in  fact,  independent. 

The  make-up  of  a  deck  for  a  MONITOR  run  is  as  follows: 

*DATE  card  (if  desired) 

*ID  card 
*XEQ  card 

Other  MONITOR  control  and  comment  cards  as  desired 

Binary  deck 

*DATA  card 

LABEL  card 

MPROB  card 

a-MATRIX  cards  (essential  for  first  problem  -  require  MALPHA  =1) 
a-VECTOR  cards  (as  necessary  when  MA>0) 

THERMO  cards 

Other  general  DATA  cards  as  needed.  Order  of  type  of  card  is  not 
important;  order  within  type  may  ba. 

Blank  card 

Subsequent  problems  starting  with  MPROB  card  and  ending  with  a 
blank  card 

End  of  file  card  (7  and  8  punch  in  column  1) 

^Indicates  MONITOR  control  cards. 
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APPENDIX  A 


EXAMPLES  OF  APPLICATIONS  OF  THE  EQUILIBRIUM  CODE 
EXAMPLE  1.  A  System  of  Gaseous  Molecules 

Consider  the  gaseous  system  composed  of  atomic  and  molecular 
oxygen  and  the  oxides  of  carbon,  with  a  known  initial  composition. 

Let  the  second  column  in  the  table  below  contain  the  list  of  species. 
The  column  headings  of  the  third  and  fourth  columns  are  the  chemical 
elements,  and  the  entries  in  these  columns  are  the  a-matrix. 


X 

o 

•^elements 

species^^^^^ 

c 

0 

al 

a2 

0 

0 

0 

1 

WBM 

1* 

0.50 

°2 

0 

2 

3 

0.50 

CO 

1 

2* 

0 

C°2 

2 

B 

4 

Q 

0.5 

1.5 

^"denotes  species  chosen  as  independent 

In  this  system,  c_  =  c_'  =  2,  and  £  =  4.  The  initial  composition  of 
the  system  is  indicated  by  the  first  column,  labelled  xq,  and  the 
elemental  composition  is  indicated  by  the  final  row  labelled  Q. 

The  selection  of  independent  species,  from  which  the  set  of 
independent  stoichiometric  reactions  follows,  can  be  made  in  any  of 
several  (in  this  case,  five)  ways.  The  obvious  restrictions  are  that 
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at  least  one  species  containing  each  element  must  be  included,  and 
that  both  0  and  0^  may  not  be  included,  since  these  two  species  have 
rows  in  the  a-matrix  which  are  multiples  of  one  another.  (In  this 
example,  these  conditions  are  redundant;  in  general  they  need  not  be.) 
The  selection  made  by  the  a-vector  in  column  five,  labelled  al, 
specifies  02  and  C02  as  independent.  The  derived  v-matrix  and  nor¬ 
malized  q-vector  for  this  a-vector  are  shown  below. 


This  formulation  of  the  equilibria  might  be  appropriate  for  temper¬ 
atures  near  2000°K,  at  atmospheric  pressure,  where  most  of  the  carbon 
in  the  system  would  be  C02  and  very  little  0  atom  would  be  present  at 
equilibrium.  The  second  a-vector,  a2,  might  be  more  appropriate  at 
much  higher  temperatures. 

EXAMPLE  2.  A  System  Containing  One  Solid  Species 

Consider  the  carbon-oxygen  system  of  Example  1  with  an  elemental 
composition  in  which  carbon  is  in  excess.  In  addition  to  the  four 
oxygen-containing  species,  let  us  also  include  gaseous  monatomic 
carbon,  C(gas),  and  solid  carbon,  C(s),  and  consider  the  possibility 
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that  C(s)  may  or  may  not  all  evaporate  at  equilibrium.  Let  the 
a-matrix,  a-vectors,  and  composition  be  given  in  a  table  like  that 
used  in  Example  1. 


"denotes  species  chosen  as  independent 
'choice  dictated  by  code 

In  order  that  its  free  energy  function  be  handled  properly  (by 
equation  22),  the  solid  species,  C(s),  must  be  listed  last  in  the 
a-matrix.  No  initial  mole  fractions  have  been  shown,  since  the  Q- 
vector  is  sufficient  for  equilibrium  computations,  and  since  it  is 
difficult  to  imagine  a  system  of  this  atomic  composition  being 
gaseous  at  ambient  temperatures. 

For  the  code  to  determine  whether  the  system  consists  of  one  or 
two  phases  at  equilibrium,  it  must  be  supplied  specifications  for  both 
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possibilities.  The  control  constants  to  be  supplied  are  0,  which 

must  be  1,  and  £_* ,  which  may  be  zero  or  one.  Two  specifications  of 

the  independent  components  are  needed.  The  specification  for  the 

two-phase  solution  is  required  by  the  code  to  have  the  solid  species 

as  the  cth  independent  component;  the  a-vector  labelled  al  is  used 

for  this  case,  and  hence  has  alg  =  c_  =  2.  The  specification  for  the 

one-phase  solution  must  have  the  solid  species  as  the  last  dependent 

species;  the  a-vector  a2  is  used  for  this  case,  and  hence  has  a2g 

=  s  =  6.  There  are  no  other  required  relationships  between  the  pair 

of  a-vectors;  the  choices  of  CO  as  an  independent  component  in  both 

cases,  and  of  C(gas)  as  the  second  component  in  a2,  are  based  on  the 

desire  to  use  species  which  will  be  present  in  appreciable  mole 

fraction  at  equilibrium  in  this  oxygen  deficient  system. 

The  two-phase  trial  solution  is  obtained  without  regard  to  the 

value  of  n  and  it  is  accepted  or  rejected  depending  on  whether 

solid 

n  . . ,  is  found  to  be  positive  or  negative.  The  one-phase  trial 
solid 

solution  is  obtained  by  treating  the  solid  as  a  gas  having  an  artifi¬ 
cially  high  value  of  F*,  thereby  limiting  the  system  to  the  first 

s 

(s  -  1)  species.  It  is  accepted  or  rejected  depending  on  whether  or 
not  the  mole  fractions  obtained  are  compatible  with  the  two-phase 
equilibria.  If  the  trial  solution  is  accepted,  the  calculation  is 
finished;  if  not,  the  other  case  is  solved  and  accepted  automatically. 
Numerical  difficulty  can  be  encountered  in  cases  where  nsolid  is 
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large  and  negative  or  in  cases  which  include  gaseous  species  having 
very  small  mole  fractions. 

EXAMPLE  3,  A  System  with  £<£* 

Let  the  system  consist  of  the  s_  =  5  gaseous  species  H,  HCN, 
CN,  and  c_’  =  3  elements  are  involved,  and  the  a-matrix  is 

then: 


However,  neither  C  nor  N  appear  in  other  than  one  to  one  combination, 
and  the  first  and  third  columns  are  redundant.  Thus  c_  =  2  must  be 
specified  in  order  for  the  code  to  handle  this  case, and  the  super¬ 
fluous  column  in  the  a-matrix  must  be  omitted. 

EXAMPLE  4.  Ionization 

Let  the  system  consist  of  the  s_  =  3  gaseous  species  He,  He+,  and 
e”.  Although  there  is  but  one  chemical  element,  in  the  ordinary 
sense,  the  desired  equilibrium.  He  ^  He+  +  e  ,  can  be  treated  by 
considering  e  to  be  an  element  whose  coefficient  in  the  formula  of 
the  positive  ion,  He+,  is  negative.  The  a-matrix  is  then  written: 
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X 

o 

He 

e 

1.0 

He 

1 

0 

0 

He* 

1 

-1 

0 

e” 

0 

1 

Q 

1.0 

0 

Let  the  initial  gas  be  pure  He,  as  indicated  by  the  column  of  initial 

mole  fractions,  x  .  Since  the  system  is  electrically  neutral,  the 
o 

resulting  value  of  Q  for  the  column  representing  e  is  zero.  This  is 

permissible  since  the  Qfs  being  zero  does  not  restrict  any  independent 

component  to  a  zero  mole  fraction. 

The  thermodynamic  functions  for  He  and  e  are  readily  handled  by 

assigning  e  a  zero  value  of  (AH,-  )  and  a  value  of  (AH^.  )  for  He 

o  o 

equal  to  the  ionization  potential  of  He.  Multiple  ionization,  nega¬ 
tive  ions,  etc.,  are  handled  by  simple  extension  of  these  principles. 
EXAMPLE  5.  Restricted  Equilibria 

4  7  8 

The  existing  general  methods  *  *  for  treating  chemical  equi¬ 
libria,  including  the  Brinkley  method  as  incorporated  in  the  HUG  code, 
are  based  on  the  specification  of  the  elemental  composition  of  each 
species.  Hence  the  system  is  constrained  only  by  the  conservation  of 
the  chemical  elements,  and  the  initial  composition  is  irrelevant.  All 
stoichiometrically  possible  reactions  are  covered  by  the  complete  set 
of  independent  equilibria  represented  by  Brinkley's  v -matrix. 
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It  is  sometimes  desired  to  consider  the  conditions  reached  by  a 
system  when  a  particular  group  of  reactions  has  reached  equilibrium 
and  other  reactions  either  have  not  occurred  or  have  occurred  to  some 
arbitrary  extent.  Provided,  of  course,  that  the  independent  reactions 
in  the  equilibrium  group  are  fewer  than  a  complete  set,  the  non¬ 
equilibrium  group  is  meaningful,  and  there  are  additional  constraints 
on  the  system  which  depend  on  the  initial  conditions  and/or  the  user’s 
whim.  To  consider  such  a  problem  by  the  general  equilibrium  codes,  it 
is  necessary  to  introduce  the  additional  constraints  by  altering  the 
a-matrix  and  Q-vector. 

In  many  cases  this  can  be  done  by  inspection.  Thus  it  is  trivial 
to  simply  omit  from  consideration  a  species  which  is  not  present 
initially  and  not  formed  by  any  reaction.  It  is  nearly  as  simple  to 
retain  a  species  in  the  system  but  render  it  inert  by  adding  to  the 
a-matrix  a  column  which  contains  zero  in  all  rows  except  that  for  the 
particular  species.  This  might  be  useful,  for  example,  in  considering 
equilibria  among  oxides  of  nitrogen  in  the  presence  of  diluent 
which  is  to  be  treated  as  inert. 

It  is  perfectly  possible  to  consider  systems  in  which  each 
species  is  involved  in  at  least  one  equilibrium  reaction,  but  all 
possible  equilibria  are  not  allowed.  A  simple  example  is  that  of 
allowing  the  equilibrium  0^  ^  0^  +  0  while  excluding  the  reaction 
02  20.  The  required  a-matrix  must  have  c_  =  2.  The  matrices  given 
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below  illustrate  two  of  the  many  possible  ways  of  accomplishing  this. 


The  first  of  these  retains  the  proper  ct-matrix  and  simply  augments  it; 
the  second  does  not,  but  uses  instead  a  matrix  borrowed  from  the  anal¬ 
ogous  case  of  N^O  +  0. 

One  case  of  restricted  equilibrium  which  is  frequently  considered 
in  combustion  and  detonation  problems  is  that  in  which  the  pressure- 
dependent  equilibria  are  established,  but  dissociation-recombination 
equilibria  are  not.  This,  too,  is  a  case  where  one  additional  con¬ 
straint  is  needed,  and  it  can  be  introduced  in  any  system  by  augment¬ 
ing  the  proper  ct-matrix  by  a  column  containing  all  ones.  (Any  constant 
would  do,  as  would  any  of  the  sets  of  coefficients  generated  by  the 
methods  indicated  by  Kaskan  and  Schott.  )  When  the  column  of  ones  is 
used,  it  has  a  Q  value  equal  to  ngas* 

It  is  possible  to  construct  augmented  ot-matrices  for  considering 
restricted  equilibria  achieved  by  an  arbitrary  mechanism  of  stoichio- 
metrically  bona  fide  reactions  in  an  arbitrary  system  of  species.  A 
recipe  for  doing  this  is  presented  here,  not  so  much  to  encourage  its 
use  in  preference  to  construction  by  inspection,  but  to  demonstrate 
its  existence  in  the  general  case. 
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Let  there  be  given  the  list  of  £  species,  together  with  the 
elemental  composition  (row  of  the  proper  a-matrix)  of  each,  and  the 
list  of  R  chemical  reactions  to  be  considered,  together  with  the 
stoichiometric  coefficients,  Bkt»  specifying  the  involvement  of  the 
t**1  species  in  the  k**1  reaction.  Reactants  and  products  are  dis¬ 
tinguished  by  sign. 

Let  us  use  the  symbol  [  ]  to  denote  matrices,  and  [  ](m)  to  indi¬ 
cate  their  dimensions,  n_  rows  and  in  columns.  The  proper  a-matrix  we 

s 

may  denote  by  [ac,](c’),  whose  rank  is  presumed  to  be  equal  to  c_* . 

Let  the  coefficients  gk  form  the  matrix  [g  ](s),  and  let  its  rank  be 
r_  <  R.  Extract  r_  independent  rows  from  [/3^],  and  let  these  be  denoted 
by  [g](s).  must  be  less  than  r_’  =  (s_  -  c.’ )  so  that  restrictions 
actually  exist;  the  number  of  independent  restrictions  is  (r'  -  r), 
and  this  is  the  number  of  extra  columns  to  be  constructed. 

Now  fix  the  order  of  species  in  the  following  way.  Select  for 
the  first  c_'  species  (t  *  1,  ...  cj)  a  set  whose  rows  of  [ac,]  form  a 

-j  * 

square  matrix,  [aJ  ],  whose  determinant  is  not  zero.  Let  the  corre¬ 
sponding  columns  of  [g]  be  denoted  by  [g^,],  and  the  remaining  r_' 
columns  be  [g^,].  Select  from  [g^,]  a  group  of  r  species  whose 
columns  form  a  square  matrix  [g^]  whose  determinant  is  not  zero,  and 
list  these  species  last  (jt  =  (£-r+l),  ...  s_).  The  remainder  of  [g] 
is  denoted  by  [g.]. 

That  this  can  always  be  done  is  seen  by  observing  that  a  reaction 
coefficient  matrix  for  a  complete  set  of  r'  independent  stoichiometric 
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r  * 

reactions  is  obtained  by  augmenting  Brinkley's  matrix  [v](c')  by  the 
negative  unit  matrix  of  dimension  £*  to  obtain  [v'](s  ),  and  that 
-[0^,]  is  the  matrix  by  which  one  multiplies  [v']  on  the  left  to 
obtain  [0],  Since  the  rank  of  a  product  cannot  exceed  that  of  either 
factor,  the  rank  of  [0^t]  is  not  less  than  r. 

Now  the  job  is  to  construct  (r*  -  r)  additional  columns  to 
augment  [a  ,]  to  form  the  final  matrix,  [a],  having  £  =  (£  -  r)  inde- 
pendent  columns.  Subject  only  to  the  requirement  of  independence, 
the  first  c  rows  of  [a]  can  be  assigned  arbitrarily.  Thus  it  is  con¬ 
venient  to  make  the  first  £*  rows  of  the  new  columns  zero,  and  the 
next  (c  -  c')  rows  a  diagonal  (conveniently  unit)  matrix.  The  remain¬ 
ing  £  rows  (£  =  (£+1),  . . .  £)  are  then  determined  by  the  relationship 

[a1]  =  -  [g.rte.]  Caj], 

where  [0.]”^  is  the  inverse  matrix  of  [0.],  [a^]  is  the  first  £  rows 
of  [a],  and  [a*]  is  the -last  r  rows  of  [a].  This  relationship  follows 
from  the  requirement  that  each  element  be  conserved  by  each  reaction, 
viz.  [0]  [a]  =  0,  and  although  it  applies  to  all  columns  of  [a],  it 
need  be  applied  only  to  the  newly  added  ones.  The  operation  involving 
is  of  course  simplified  by  operating  by  rows  on  the  original 
[0]  until  [0.]  is  the  unit  matrix  of  dimension  r.  In  order  to  ac¬ 
commodate  the  loading  format  for  the  a-matrix,  the  added  columns  must 
be  cleared  of  fractions. 
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If  the  computation  is  to  be  restricted  to  states  attainable  from 
the  initial  composition  by  the  equilibrium  reactions  only,  the  proper 
values  of  Q  are  generated  by  equation  (16),  whether  the  columns  of  [a] 
represent  true  chemical  elements  or  pseudo-elements.  An  arbitrary 
nonzero  extent  of  any  nonequilibrium  reactions  may  also  be  considered 
by  using  an  improper  set  of  Q's  derived  via  equation  (16)  from  a  set 
of  nt's  (see  XMOLE  CARDS)  which  represent  the  formal  composition 
obtained  by  adjusting  the  initial  composition  for  the  arbitrary  re¬ 
actions.  Such  a  set  of  Q's,  which  may  be  inconsistent  with  that  based 
on  XCOMPO  data,  may  be  used  by  loading  it  via  a  QATOM  card.  (The 
inconsistency  arises  from  the  restriction(s)  on  the  allowed  equilibria 
and  does  not  represent  a  violation  of  the  elemental  mass  conservation 
requirements. ) 
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APPENDIX  B 


A  LIBRARY  OF  THERMODYNAMIC  COEFFICIENTS  DERIVED  FROM  THE  JANAF  TABLES 
Prepared  by  Charles  W.  Hamilton 

Tables  B.1-B.3  list  a  library  of  coefficients  used  for  calcu¬ 
lating  thermodynamic  functions  of  a  number  of  species  by  equations 
(17)-(20).  The  unit  of  enthalpy  in  equations  (17)  and  (18)  is 
calories  per  gram  mole;  the  left  hand  sides  of  equations  (19)  and  (20) 
are  dimensionless.  The  unit  of  temperature  is  °K  throughout. 

The  JANAF  tables^  were  sole  reference  for  the  enthalpies,  free 
energies,  heats  of  formation,  and  specific  heats  used  in  generating 
these  coefficients. 

Table  B.l  contains  the  values  of  a,  b,  c,  and  d  in  equation  (18) 
for  the  temperature  range  200-400°K  found  by  simultaneous  solution  of 
the  four  equations: 

a  +  bT  +  cT2  +  dT3  =  H'(T),  T  =  200,  300,  400°K 

b  +  2cT  +  3dT2  =  Cp(T),  T  =  300°K 

where  H'(T)  =  H(T)  -  H(0°K),  and  Cp  is  the  tabulated  specific  heat. 

Tables  B.2  and  B.3  contain  values  of  the  coefficients  in 
equations  (19)  and  (20)  for  the  two  overalpping  ranges  of  elevated 
temperatures,  600-3000°K  and  1600-6000°K,  respectively.  The  values  of 
A,  B,  C,  D,  and  E  were  found  by  least  squares  fits  of  equation  (19)  at 
100°  intervals.  Equation  (20)  was  then  solved  for  k  at  each  temper- 
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ature,  and  the  arithmetic  mean  value  of  the  k's  thus  generated  was 
tabulated. 

Heats  of  formation  at  0°K  taken  directly  from  the  JANAF  tables 
are  also  listed  in  tables  B.2  and  B.3. 


63 


TABLE  B.l 


TEMPERATURE 


SPECIES 

RANGE 

d 

c 

BF3 

200-400 

-6.600-06 

♦1.538-02 

+4.673+00 

♦2.008+02 

B2H6 

200- MOO 

♦7.000-06 

♦1.225-02 

♦4.135+00 

♦3.240+02 

HER 

200-400 

+0.000+00 

♦5.000-05 

+6.935+C0 

-4.900+00 

BR2 

200-400 

-2.900-06 

+3.900-03 

♦7.090+00 

-5.690+01 

HCN 

200-400 

-3.500-06 

+7.890-03 

♦4.849+00 

+1.574+02 

CH4 

200-400 

+1.050-05 

-5.140-03 

♦8.790+00 

-4.490+01 

CO 

200-400 

♦  1.500-06 

-1.340-03 

♦7.370+00 

-4.390+01 

C02 

200-400 

-7.000-06 

+1.234-02 

+3.409+00 

+3.124+02 

C2H2 

200-400 

-1.050-05 

♦1.839-02 

♦2.399+00 

♦3.244+02 

C2H4 

200-400 

♦7.400-06 

♦4.040-03 

+5.554+00 

♦2.434+02 

C2H40 

200-400 

♦1.050-05 

♦6.500-03 

♦4.775+00 

♦3. 160+02 

C2N2 

200-400 

-1.190-05 

♦1.830-02 

♦5.850+00 

-3.490+01 

C302 

200-400 

-1.640-05 

♦2.715-02 

+3.910+00 

♦1.000+02 

C4N2 

200-400 

-2.  100-05 

+3.224-02 

+6.929+00 

-1.335+02 

HCL 

200-400 

+1.000-07 

-3.000-05 

♦6.961+00 

-9.300+00 

CL2 

200-400 

-3.300-06 

♦5.210-03 

+5.911+00 

+5.860+01 

HF 

200-400 

-3.000-07 

+3.600-C4 

♦6.856+00 

-1.030+01 

F2 

200-400 

-2.000-07 

♦2.320-03 

+6. 187+00 

+6.720+01 

SF6 

200-400 

-3.950-05 

♦6.404-02 

-4.365+00 

♦  7. 144+02 

HI 

200-400 

+1.000-07 

♦1.00C-05 

♦6. 936+00 

-2.300+00 

H2 

200-400 

-4.800-06 

♦5.360-03 

♦5.001+00 

♦1.866+02 

H20 

200-400 

♦  1.800-06 

-1.060-03 

♦8. 183+00 

-2. $10+01 

H202 

200-400 

-3.640-05 

+4.290-02 

-5.554+00 

+1.404+03 

H2S 

200-400 

+2.400-06 

-7.500-04 

♦7.984+00 

+4.400+00 

NH3 

200-400 

♦4.100-Q6 

-1.230-03 

+8.021+00 

♦4.600+00 

12 

200-400 

-1.600-06 

+2.230-03 

+7.938+00 

-1.021+02 

NO 

200-400 

♦2.800-06 

-2.710-03 

♦8.008+00 

-2.310+01 

N02 

200-400 

♦1.800-06 

+2.130-03 

♦  7.  123+00 

♦7.580+01 

N2 

200-400 

♦9. 000-07 

-8.000-04 

♦7.204+00 

-2.750+01 

N20 

200-400 

-4.600-06 

♦9.630-03 

♦4.738+00 

♦1.458+02 

N203 

200-400 

-9. 100-Q6 

♦1.783-02 

+8.493+00 

-5.910+01 

N204 

200-400 

-1.320-05 

♦2.712-02 

+5. 767+00 

+  1. 252  +  02 

N205 

200-400 

-1.520-05 

♦3.372-0 2 

♦4.677+00 

+1.792+02 

02 

200-400 

+7.000-07 

♦1.700-04 

+6.732+00 

♦3.420+01 

S02 

200-400 

-1.000-07 

♦4.530-03 

♦6.883+00 

♦7.280+01 

03 

200-400 

+5.000-07 

♦5.000-03 

+6.265+00 

♦1.480+02 

S03 

200-400 

-  1.790-05 

♦2.760-02 

♦5.400-01 

+6.540+02 
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TABLE  B.2  (continu'd) 
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TABLE  B.3  (oontinuad) 
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